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FOREWORD 


The  program  described  in  this  report  reflects  the  research 
phase  of  the  numerical  facilitation  of  the  hybrid  stress  finite 
element  method  for  the  large  deflection  stress  analysis  of 
multi-layered  aircraft  tires.  The  work  was  administered  by  the 
Air  Force  Flight  Dynamics  Laboratory,  WPAFB,  Ohio,  under 
Contract  Nos.  F33615-72-C-1693  and  F33615-73-C-3002  for  the 
period  of  10  January  to  12  November  1973  under  Project  1369, 
"Mechanical  Sub-Systems  for  Advanced  Military  Flight  Vehicles," 

Task  No.  136903,  "Landing  Gear  System  Ground  Contact  Components 
for  Advanced  Military  Flight  Vehicles."  Dr.  H.  K.  Brewer 
served  as  the  principal  technical  monitor  for  the  Air  Force. 

The  authors  are  indebted  to  Marianne  M.  Montgomery  whose 
insight  into  the  problems  of  large-scale  computer  program 
development  has  allowed  the  completion  of  the  work  leading 
up  to  this  program. 

The  contractor  report  number  is  MSNW  73-303-1. 

This  technical  report  has  been  reviewed  and  is  approved. 

j?r,  2£  

Kennerly  H.  Digge^k) 

Chief,  Mechanical  Branch 

Vehicle  Equipment  Division 

Air  Force  Flight  Dynamics  Laboratory 


ABSTRACT 


Presented  is  a description  of  the  FORTRAN/COMPASS  computer  code 
for  the  large  deflection  stress  analysis  of  multi-layered  aircraft 
tires.  The  program  is  modulated  into  nine  overlays  within  the  frame- 
work of  dynamic  storage  allocation  and  is  operational  on  the  CDC-6600 
machine  under  the  SCOPE  3.3  system. 
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1.  INTRODUCTION 


The  computer  code  for  the  stress  analysis  of  aircraft  tires 
is  designed  to  solve  the  following  problems: 

• Inflation  of  the  lifted  but  unloaded  tire 

• Rotation  of  an  inflated  but  unloaded  tire 

• Contact  problem  of  a statically  loaded  and  inflated  tire 

• Contact  problem  of  a loaded  rolling  aircraft  tire. 

The  code  is  subdivided  into  eight  overlays  within  the  framework 
of  dynamic  storage  allocation.  In  the  first  four  overlays  the  input 
data  is  reduced  to  set  up  quantities  associated  with  the  geometrical 
configuration.  The  fourth  overlay  calculates  the  element  stiffness 
and  load  matrices,  which  are  assembled  in  the  fifth  overlay.  The  sixth 
overlay  contains  a direct  equation  solver  with  one  right-hand  side.  For 
the  contact  problem,  the  seventh  overlay  generates  the  flexibility  matrix 
coefficients  using  a direct  multiple  right-hand  side  equation  solver. 

The  actual  contact  problem  algorithm  is  contained  in  the  eighth  overlay. 

In  the  following  section,  the  structure  and  modulation  of  the 
code  will  be  described  in  detail. 

1,1.  Storage  Allocation  and  Input/Output  Characteristics 

The  data  management  of  the  computer  code  incorporates  those  pri- 
mary features  of  the  CDC  6600  system  which  are  necessary  for  the  efficient 
flow  of  large  sets  of  information.  Information  storage  and  retrieval 
procedures  were  designed  to  minimize: 
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• Central  memory  required 

• Input/output  access  time 

• Program  maintenance  and  modifications. 

The  following  main  features  of  the  CDC  6600  system  were  used 
to  achieve  the  above  objectives: 

• Random  access  input/output  subroutines 

• Unblocked,  unbuffered  files 

• Blank  common. 

The  random  access  subroutines  are  library  input/output  routines, 
supported  by  CDC,  which  provide  the  capability  for  direct  storage  and 
retrieval  of  records  on  a file,  as  opposed  to  sequential  files  which 
require  accessing  the  records  preceding  the  desired  one.  In  the  computer 
code,  these  mass  storage  routines  are  used  extensively  to  store  the 

• Input  data  and 

• Computed  data 

between  overlays,  which  allows  the  program  to 

• Select, 

• Input  and 

• Output 

only  those  data  which  are  necessary  in  the  particular  overlay  under 
consideration. 

Unblocked,  unbuffered  files  are  used  to  store  the  intermediate 
data;  for  instance,  in  the  calculation  of  the  element  stiffness  and  ele- 
ment load  matrix.  These  files  are  accessed  repeatedly  in  a sequential 
manner.  These  unblocked,  unbuffered  files  are  efficient  for  reading  and 
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writing  large  records,  since  the  information  is  read  directly  from  the 
disk  into  the  program  array  area.  Note  that  the  records  in  blocked  and 
buffered  files  would  first  be  read  into  an  intermediate  system  storage 
of  the  central  memory  area,  and  then  transferred  into  the  program  array 
area. 

The  proper  use  of  blank  common  allows  the  code  to  have  a general 
work  area  available,  whose  length  depends  only  on  the  field  length  de- 
clared on  the  job  card.  This  area  is  dynamically  divided  among  the 
arrays  needed  in  executing  the  current  overlay. 

Thus,  the  central  memory  disposition  of  the  codes  have  the  fol- 
lowing structure: 
word  1 


field  length) 

The  above  construction  allows  the  user  to  specify  a field  length 
tailored  to  the  data  size. 

In  the  current  code,  each  overlay  determines  the  length  of  the 
arrays  used  and  stacks  them  nose— to-tail  in  blank  common. 

The  current  code  uses  no  tapes.  It  is  realized  that  this  feature 
is  essential  in  modern  technology  and  thus  we  propose  to  perform  all 
improvements  in  the  "no  tape"  philosophy. 


1.2.  Fortran  Extended  Code 


Two  versions  of  the  code  are  provided,  one  produced  via  the  RUN 
compiler  and  the  other  the  fortran  extended  (FTN)  version.  The  RUN 
version  varies  from  the  FTN  version  by  its  IF  UNIT  tests,  presence  of 
RETURN  statements  in  overlay  main  programs  and  its  compass  decks. 

All  of  these  are  provided. 

1.3.  Library  Routines 

Besides  the  standard  library  routines,  the  codes  employ  the 
following  special  features  of  the  CDC  6600  library: 

• BUFFER  IN 

• BUFFER  OUT 

• READMS 

• WRITMS 

1.4.  Assembly  Language  Subroutines 

The  CDC  6600  assembly  language,  COMPASS  (COMPrehensive  ASSembly 
language),  is  particularly  suited  to  substanially  reducing  the  compu- 
tation time  of  looped  operations.  The  improvements  are  realized  by: 

• More  efficient  retrieval  of  array  elements 

• Overlapping  of  data  storage  and  retrieval  from  central  mem- 
ory with  multiplication,  addition  and  subtraction 

• Efficient  use  of  the  instruction  stack  which  holds  seven 
words  (up  to  28  instructions)  in  the  central  processor. 

Well  coded  compass  routines  will  execute  computational  do  loops 
from  5 to  6 times  faster  than  normal  Fortran  IV  on  the  CDC  6600. 
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In  the  present  code,  six  of  the  heavily  used  matrix  manipulation 
subroutines  are  written  in  COMPASS! 

• MATMPY 

• MATADD 

• MATSMP 

• INPRDS 

• VECMAT 

• EMULT 

'j’heje  are  six  special  compass  subroutines  to  perform  tasks  which 
standard  Fortran  IV  is  not  designed  to  handle.  Th.*.se  are  listed  below: 

• KFL 

• SSZERtj) 

■ MSTG 

• GET 

• PUT 

• STRM<)>V 

KFL  is  a COMPASS  subroutine  which  retrieves  the  field  length 
requested  by  the  job.  This  information  allows  the  program  to  use  all 
central  memory  available.  Furthermore,  for  each  data  set,  a minimum 
field  length  requirement  may  be  tailored. 

SSZERtf)  is  designed  to  set  array  values  to  zero  during  execution 
in  a minimum  amount  of  time.  The  routine  is  used  extensively  throughout 
the  code. 

MSTG,  GET,  PUT,  and  STRM(f>V  are  COMPASS  routines  which  perform 
character  and  string  manipulation.  They  are  used  in  dynamic  storage 
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allocation  and  by  the  free  field  input  reader  which  rec1  ices  considerably 
the  time  requred  by  the  user  to  enter  and  debug  his  input  data. 


2.  PROGRAM  ORGANIZATION 


The  program  is  organized  according  to  the  following  problem 

types : 

• inflation 

• inflation  and  rotation 

• inflation  and  contact 

• inflation,  rotation  and  contact 

The  general  data  flow  is  shown  below,  where  the  numbers  refer 
to  the  overlays  under  consideration. 


INFLATION 


ROTATION 


CONTACT 


In  the  above  data  flow  the  symboll  [refers  to  the  incremental 
inflation  process.  The  symbol  /^refers  to  the  contact  iteration  al- 
gorithm. 


3.  OVERLAY  DESCRIPTION 


3.1.  Overlay  (KTIRE , 0,0) 

This  overlay  controls  the  general  data  flow  as  described  in 
Section  2.  It  performs  the  initialization  of  labeled  common  blocks, 
opens  random  access  mass  storage  files  and  facilitates  the  storage 
requirements  within  the  framework  of  dynamic  storage  allocation.  This 
overlay  also  contains  various  utility  programs  and  assembly  language 
routines  for  vector  and  character  manipulations. 

3.2.  Overlay  (KTIRE,  1,0) 

All  the  input  data  are  read  from  cards  in  this  overlay  and  then 
they  are  placed  on  random  access  mass  storage  files.  The  input  data 
are  checked  for  logical  errors  which  are  summarized  at  the  end  of  the 
data  processing  phase,  using  the  subprograms  RANGE,  WRDCHK  and  COMPCHK. 
The  primary  control  parameters  are  also  set  up  here  in  the  labelled  com- 
mon blocks  SIZE  and  CONTACT. 

3.3.  Overlay  (KTIRE,  2,0) 

This  overlay  performs  the  preliminary  nodal  calculations  such  as 

• cartesian  and  curvlinear  coordinates 

• surface  vectors  of  the  undeformed  reference  surface 

• cord  angle  distribution  along  the  meridian. 

The  resulting  data  are  placed  on  random  access  mass  storage  files 

3.4.  Overlay  (KTIRE,  3,0) 

The  intrinsic  reference  element  properties,  such  as 
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element  area 


• 

o element  centroid 

• local  unit  vectors 

• local  element  vertex  coordinates 

• average  cord  angles 

are  generated  here  and  then  placed  on  random  access  mass  storage  files. 
3.5.  Overlay  (KTIRE , 4,0) 

This  overlay  sets  up  the  element  stiffness  matrix  and  load 
vector  according  to  the  hybrid  stress  finite  element  formulation  out- 
lined in  [1]. 

In  principle,  the  element  complementary  energy  matrix  and  the 
boundary  work  by  the  stress  resultants  are  calculated  here,  followed  by 
the  elimination  of  the  undetermined  stress  coordinates  using  an  out-of- 
core Choleski  decomposition  algorithm. 

In  particular,  it  performs  the  calculation  of  the 

• element  flexibility  matrix 

• element  flexibility  vector 

• homogeneous  incremental  flexibility  matrix 

• incremental  element  flexibility  vector 

• particular  incremental  element  flexibility  matrix 

• hybrid  element  load  matrix 

• hybrid  element  load  vector 

• incremental  hybrid  element  lead  matrix. 
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Since  all  these  element-wise  calculations  involve  out-of-core 


processing,  the  theme  of  fourth  overlay  is  established  by  the  Choleski 
inversion  routine  and  best  explained  by  considering  the  main  calling  se- 


quences as  follows. 

In  Table  1 the  key  subroutines  have  the  functions 

• CUBRE  sets  up  the  Gaussian  weights/nodes 

• TRCALC  sets  up  the  transformation  matrix  in  the  local 

coordinate  system  for  the  lamina  constitutive  rela- 
tions 


• DCALC 
SCALC 


calculate  the  lamina  compliance  in  a principal 
frame  from  the  constitutients  elastic  properties. 

• ENER  performs  the  calculation  of  the  complementary  energy 

matrices  as  indicated  by  Table  1. 

• HBMERGE  merges  the  layer  complementary  energy  matrices  due  to 

homogeneous  stress  field  to  obtain  the  element  flexi- 
bility matrix  Cao,  as  described  in  Section  3.1.10  [1]. 
pp 

• HAMERGE  merges  the  layer  complementary  energy  matrix  due  to 

homogeneous  and  particular  stress  field  to  obtain  the 
corresponding  element  complementary  energy  matrix  C^. 

• PQMERGE  merges  the  incremental  layer  complementary  energy 

matrices  to  form  the  corresponding  element  complementary 

energy  matrices  ACfi  and  AC 

Pq  ^q 

• WORK  calculates  the  work  done  by  the  stress  resultants  on  the 

reference  surface  displacements. 

• PRESF  sets  up  the  external  load  vectors 
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• CHOLESK 
DECOMP 
FORSUB 
BAKSUB 


performs  the  inversion  of  the  element  flexibility 

matrix  COD  using  a Choleski  decomposition  procedure 
pp 

which  takes  into  account  the  banded  symmetric  and 
staircase  structure  of  this  matrix. 


• ELMAT  forms  the  element  stiffness  matrix  and  load  vector 

according  to  Equations  (3.1.16)  of  [1]. 

• CAHSTR  perform  the  area  integrations  using  numerical 

I cubatures  as  described  in  Section  3,  Vol.  I. 

DSBQ 


Note  that  a single  pass  at  the  fourth  overlay  sets  up  the  element 
stiffness  and  load  matrices.  After  assembling  these  in  the  augmented 
structure  stiffness  matrix  and  solving  for  the  general  displacements,  one 
must  sweep  overlay  four  the  second  time  to  fetch  the  stress  coordinates, 
using  Equation  (3. 1.1. 3)  of  [1]  to  obtain  the  layerwise  stress  distribu- 
tion for  each  element. 

3*6.  Overlay  (KTIRE,  5,0) 

This  overlay  constructs  the  structure  stiffness  matrix  for  all  the 
problems  under  consideration,  and  the  structure  load  vectors  for  inflation 
and  rotation  problems. 

The  corresponding  merge  routine  takes  into  account  the  banded  and 
symmetric  properties  of  the  structure  stiffness  matrix  which  is  augmented 
by  a single  load  vector.  According  to  the  problem-size  under  consideration, 
this  augmented  matrix  is  subdivided  into  row-wise  blocks  for  out-of-core 
processing.  Each  block  is  then  placed  on  a random  access  file,  so  that 
the  element  stiffness  matrix  and  load  vector  components  may  effectively  be 
placed  in  the  appropriate  block  under  consideration. 
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After  assembling,  the  homogeneous  displacement  boundary  conditions 
are  applied  by  zeroing  out  the  corresponding  rows  and  columns  of  the 

augmented  structure  stiffness  matrix  and  placing  a finite  number  in  the 
diagonal . 

For  the  contact  problem,  the  same  merge  procedure  is  used,  however, 
the  structure  stiffness  matrix  is  augmented  by  multiple  right-hand  sides 
to  obtain  the  appropriate  flexibility  matrix  components  required  for  the 
facilitation  of  the  contact  analysis. 

3.7.  Overlay  (KTIRE,  6,0) 

A standard  Gaussian  elimination  routine  is  contained  here  which 
takes  into  account  the  banded  and  symmetric  properties  of  the  structure 
stiffness  matrix,  which  is  set  up  according  to  Section  3.6.  The  appro- 
priate inner-product  operations  are  coded  in  COMPASS  as  exhibited  by  the 
EMULT  subroutine. 

This  overlay  is  assessed  at  each  incremental  step  during  inflation 
and  also  for  the  rotation  problem.  After  each  pass,  the  resulting  solu- 
tion vector  is  used  to  update  the  geometrical  configuration  according  to 
the  intrinsic  initial  stress  formulation. 

3.8.  Overlay  (KTIRE,  7,0) 

This  overlay  sets  up  the  appropriate  flexibility  matrix  coefficients 
for  the  contact  problem,  which  are  obtained  by  a direct  multiple  right-hand 
side  equation  solver,  called  SOLVMOR.  The  relevant  inner  product  opera- 
tions in  this  routine  are  coded  in  assembly  language  with  double  overlapping, 
EMULT. 
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For  the  contact  problem  the  global  coordinate  system  is  a 
cartesian  reference  frame.  The  appropriate  element  stiffness  and  load 
matrices  are  calculated  in  the  fourth  overlay.  The  structure  stiffness 
matrix,  as  for  the  inflation  or  rotation  problem,  is  again  constructed 
in  the  fifth  overlay. 

The  seventh  overlay  reads  the  above  structure  stiffness  matrix 
and  augments  it  with  appropriate  unit  vectors  defined  by  the  candidate 
contact  nodes  under  consideration.  The  corresponding  flexibility  matrix 
coefficients  are  then  obtained  from  SOLVMOR,  which  destroys  the  structure 
stiffness  matrix  during  the  reduction  process.  Thus,  during  the  contact 
iteration,  the  structure  stiffness  matrix  is  re-fetched  from  the  mass 
storage  files  created  by  the  fifth  overlay. 

3.9.  Overlay  (KTIRE,  9,0) 

This  overlay  performs  the  iterative  contact  analysis  as  described 
in  Section  3.2  [1] . 


4.  FILES  AND  COMMON  BLOCKS 


The  program  control  variables  are  transmitted  via  labelled  com- 
mon blocks,  while  the  relevant  input  or  calculated  data  are  placed  on 
sequential  or  random  access  mass  storage  files. 

4.1.  Common  Blocks 


BCINDEX  contains  the  information  for  dynamic  storage  allocation: 

• NOPOS  or  number  of  positions  available  for  allocation. 

• KSPACE  or  space  available  for  allocation. 

• INDEX  or  pointers  to  currently  defined  arrays 


CONTACT  consists  of  the  information  for  the  iterative  contact  algorithm: 

• NORING  or  number  of  contact  rings 

• KRING  or  current  ring  number 

• NODMAX  or  maximum  number  of  nodes  in  a ring 

• KSZAX  or  total  number  of  candidate  contact  nodes 


ERROR 


accounts  for  all  logical  input  data  errors  and  includes  the 
following  informations : 

NERR  or  running  count  of  input  errors 

NERRS  or  record  numbers  of  those  with  errors 

NERLIM  or  maximum  number  of  errors  to  be  counted 

KERR  or  indicator  as  to  whether  current  block  has  errors. 


FILES  contains  the  names  of  all  data  files  defined  by  Fortran  IV 

Hollerith  form. 
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INDTA 


is  comprised  of  information  about  the  last  record  read  by  the 
input  reader: 

NWRD  or  number  of  items  present 
ITYP  or  type  of  each  item  present 
NREC  or  record  number 
NCRD  or  card  number 

DTA  or  value  of  each  item  in  the  record 


MATS1Z 


consists  of  the  structure  stiffness  matrix  characteristics: 
NUMBK  or  number  of  blocks 
NBKSI  or  block  size 

NMIQ  or  bandwidth  including  right-hand  side. 

NPB  or  nodal  points  per  blocks 

NEQ  or  number  of  equations  per  blocks 

NMAX  or  total  number  of  equations 

NORHS  or  number  of  right-hand  sides 

HRBKSI  or  the  size  of  the  flexibility  matrix  block 


PRINTS 


RECORD 


controls  all  optional  intermediate  printing: 

KPRINT  or  an  array  indicating  which  intermediate  values  the 

user  wants  printed 

LINLIM  or  the  maximum  number  of  single  spaced  lines  pe  p g 

the  variable  names  for  all  named  random  access  records. 


contains 
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RETRIV 


consists  of  information  to  determine  the  blocksize  of  the 
structure  stiffness  matrix: 

• LENCOM  or  the  address  of  the  beginning  of  blank  common  for 
each  overlay. 


SIZE 


’] 

is  comprised  of  the  input  control  parameters : 

• NOEL  or  number  of  elements 

• NNODE  or  number  of  nodes 

• RADIUS  or  radius  of  meridian  reference  curve  rotation 

• NRHO  or  number  of  fitting  coefficients  for  the  meridian 
reference  curve 

• NPRHO  or  number  of  data  p )ints  for  the  meridian  reference 
curve 

• MAXLAY  or  maximum  number  of  layers 

• GREEN  or  green  angle 

• SPEED  or  rotational  speed 

• INCR  or  number  of  increments  for  nonlinear  inflation. 
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4.2. 


Sequential  Files 


All  the  values  contained  in  the  files  described  below  are  directly 
calculated  in  the  fourth  overlay.  File  descriptions  are  given  in  Table  2 


Table  2 


Sequential  File  Description 


File  Name 

Contents 

Size 

How  Many 
Generated 

KBMAT 

B matrix  used  in 
ELMMAT 

Maximum 
315  words 

One  record  per 
layer 

KBQMAT 

BCBQ  matrix  used 
in  ELMMAT 

Maximum 
315  words 

One  record  per 
layer 

KHAMAT 

HP  matrix  used 
in  ELMMAT 

Maximum 
189  words 

One  record  per 
layer 

KHBMAT 

H matrix  to  be 
inverted  via 
CHOLESK  inversion 

Maximum 
378  words 

One  record  per 
layer 

4.3.  Random  Access  Files 

Following  is  a description  of  all  records  of  all  random  access 
files.  The  fourth  table  entry  indicates  whether  the  information  is  in- 
put or  calculated  by  the  code. 
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Table  3 


Random  Access  Files 


File  Name 

Record 

Name 

Contents 

I/C 

Size 

KCTMAT 

KCTC 

Contact  data 

I 

NORING  * 
(NODMAX  + 2) 

KCTMAT 

KCTC2 

X Coordinate  Array- 
Contact  case 

C 

KSZAX 

KCTMAT 

KCTC  3 

Contact  Case  ALPHA 
array 

C 

Maximum 

KSZAX 

KCTMAT 

KCTC  5 

Total  Nodal  Load 
Vector-Contact  Case 

C 

Maximum 

KSZAX 

KHMAT 

Numbered 

one/layer 

Decomposition  of  the 
HB  matrix 

C 

Maximum  378 

KLADAT 

Numbered 
one /element 

Element  Layer  Data 

I 

Maximum 
MAXLAY  *10 

KONDAT 

KAR 

Elements  areas 

C 

NOEL*l 

KONDAT 

KAV 

Average  Cord  Angles 

C 

NOEL*l 

KONDAT 

KB  ETA 

Curvefit  betas 

I 

NPBETA*2 

KONDAT 

KCA 

Cord  Angles 

C 

NNODE*l 

KONDAT 

KCAR 

Cartesian  coordinates 

C 

NNODE*3 

KONDAT 

KCARC 

Cartesian  Data 

I 

NNODE*6 

KONDAT 

KCEN 

Centroid 

c 

NOEL*3 

KONDAT 

KCLC 

Curvilinear  coordi- 
nates 

I 

NNODE*2 

KONDAT 

KDIS 

Displacements 

I 

NNODE*6 

KONDAT 

KELN 

Elements'  Nodes 

I 

NOEL*3 

KONDAT 

KFOR 

Forces 

I 

NNODE*3 

KONDAT 

KINCR 

Nonlinear  Increment 
Data 

I 

INCR*2 

KONDAT 

KIPV 

Inplane  Vertex 
Coordinates 

c 

NOEL*6 

KONDAT 

KLUV 

Local  Unit  Vectors 

c 

NOEL*9 

KONDAT 

KRHO 

Curvefit  Rhos 

I 

NPRHO*2 
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File  Name 


Record 

Name 


Contents 


I/C 


Size 


E 


KONDAT 

KPMAT 


KSV 


Surface  Vectors 


KRMAT 


Numbered  Upper  half  element 
one/element  load  & stiffness 
matrix 
,-l 


C 

C 


Numbered 

one/layer 


The  H * matrix  after 
forward  substitution 


KSLMAT 

KSTFIL 


KSLV 


Solution  vector 


Numbered-one 
per  merged 
block 


Output  merged  K 
matrices 


C 

C 


Numbered 
one /block 


SOLVMOR  solution 
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NN0DE*9 
135  words 


Maximum 

(NOLAY*21-7) 

*21 


NN0DE*5 


Maximum  NEQ* 
(KBAND+KSZAX) 


Maximum 
NEQ*  KSZAX 


4.4.  Data  Flow  Chart 


CARO  INPUT 


PROGRAM  FLOW  CONTROL 


KELN 

KFOR 

KOIS 

KRHO 

KCLC 

KBTA 

KINCR 

KCARC 


INPUT  & INITIALIZATION 


KONDAT, 


KCTC 

KCTC3 


KCTMAT, 


KLADAT, 


PRELIMINARY  NODAL  CALCULATIONS 


KONDAT, 


KONDAT, 


NONLINEAR  LOOP 


KONDAT, 


PRELIMINARY 


LEMENT  CALCULATIONS 


KONDAT 


GENERATES  El 


EMENT  STIFFNESS  MATRICES  AND  ELEMENT 
LOAD  MATRICES,  ONE  PER  ELEMENT 


KELN 

KSV 

KLUV 

KIPV 

KAV 

KFOR 

KCA 

KINCR 

KCEN 


KONDAT, 


Table  4 


^KLADATrJ 


1 


\ BEGIN  element  loop 


4 


V 


\ REWIND  KHAMAT 
\ REWIND  1'HBMAT 

I rewind  kbqmat 


\ BEGIN  LAYER  LOOP 
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DECOMP 


FORSUB 


BAKSUB 


KB  MAT 


KB  MAT 


END  ELEMENT  LOOP. 

I LOOP  ONCE  PER  ELEMENT 


REWIND  KPMAT 


Merge  Structure  Stiffness  and  Load  Matrices 


REWIND  KPMAT 


MERGE 


CALCULATE  SOLUTION  VECTOR 


/^KONDAT  WL 


( KCTMAT 


> 


rrcTzr^ 


./Are  ^ 
<^we  doing 


sjinear 

x 

(No) 


/i\re 

<W  doing  , 
r otHion  a1/3se  ? 


/ Call  over 
lays  3,4,5, 

& 6 for  ro- 
tational ca4e 


<we  doing 
c<*ntact  i 


Calculate  Flexibility  Matrix 


KCTC 

KCTC5 


SOLVMOR 


Brings  X Coordinates  of 
current  ring  into  con- 
tact plane  and  updates 
wheel  load. 


KCTC 

KCTC2 

KCTC3 

KCTC5 


CONTACT 


KCTMAT^ 

1 

5.  INPUT  DATA  PREPARATION 


tkgkj 


There  are  no  inherent  units  assigned  to  any  of  the  variables. 
Thus,  the  user  is  free  to  use  any  desired  set;  however,  once  defined, 
the  units  must  be  consistent  for  all  parameters. 

Format-free  processing  is  employed  for  both  floating  point  and 
integer  variables.  The  current  character  manipulation  routine  does 
not  allow  blanks.  For  instance,  the  following  floating  point  numbers 

.1567  or  -.1568 

are  not  recognized.  They  must  have  the  form 

0.1568  or  -0.1568 

Numeric  specifications  must  be  separated  from  each  other  by  at 
least  one  blank  column.  Thus  the  following  format  is  accepted: 

1.  0.5  0.2  -0.1 

5.1  Title  Cards 

No  restriction  is  imposed  on  the  number  of  title  cards.  Each 
of  these  cards  must  start  with  a slash  in  column  1.  Example: 

/ DATA  SET  NUMBER  1 
/ PREPARED  10/1/72 
/ AIRCRAFT  TIRE 
/ 
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- .Mill 


5.2  Control  Cards 


A name  is  assigned  to  each  of  the  control  variables.  The  name 
is  followed  by  the  appropriate  numerical  specification.  The  last  char- 
acter on  each  of  these  cards  must  be  the  slash  symbol,  which  is  separated 
from  the  words  or  numeric  specifications  by  any  number  of  blank  columns. 

These  control  cards  follow  the  schedule  below: 


BEGIN  DATA  INPUT  / 

BEGIN  CONTROL  PARAMETERS  / 
NODES  N / 

ELEMENTS  E / 

RADIUS  R / 

NRHO  M / 

NPRHO  N / 

LAYERS  L / 

GREEN  a / 

SPEED  0)  / 

CONTACT  AW/ 

NONLINEAR  / 

NORING  NR  / 

NODMAX  NM  / 

INCREMENTS  IN  / 

END  CONTROL  PARAMETERS  / 

/ 

/ 


The  above  numeric  control  variables  are  defined  as  follows 
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N = Number  of  nodes  (integer) 

E = Number  of  elements  (integer) 

R = Rotation  radius 

M = Number  of  curvefitting  parameters  for  the  median  section  (integer) 

N = Number  of  prescribed  data  points  for  the  meridian  section  (integer) 

L = Max  .mum  layer  number  (integer) 

a = Green  angle 

to  = Rotational  speed 

A = Initial  deflection  for  contact 

NR  = Number  of  rings 

NM  * Maximum  number  of  nodes  in  a ring 
IN  = Number  of  increments  for  nonlinear  inflation 
W = Maximum  wheel  load 
5.3  Nodal  Data 

The  node  number  and  the  corresponding  curvilinear  coordinates 

f 

are  specified  here. 

BEGIN  NODAL  DATA  / 

i 0 

END  NODAL  DATA  / 

/ 

/ 

The  numeric  variables  i,0,  and  cf)  are 
i = Node  number  (integer) 

0=  Parallel  in  radians  (floating  p.) 

<f>  = Meridian  in  radians  (floating  p.) 
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5.4  Element  Data 


There  are  three  nodes  associated  with  each  surface  element. 

The  node  assignments  must  follow  the  right-hand  rule  according  to  the 
outward  normal  direction.  Furthermore,  the  first  node  number  must  be 
the  smallest . 

Within  this  record,  the  elastic  constants  are  also  specified  for 
each  layer  within  the  element,  in  a principal  reference  frame.  Thus, 


BEGIN 
E j 

ELEMENT  DATA 
N1  N2 

/ 

N3 

N4  1 

L 1 
• 

• 

t r 

** 

VR  Ec 

Vc  b 

Sc  SR  f / 

• 

L k t r 

END  ELEMENT  DATA  / 

VR  Ec 

vc  b 

Sc  SR  £ ' 

The  numeric  characters  in  the  above  records  are 
j = Element  number  (integer) 

N^,  N^,  = Node  numbers,  assigned  according  to  the  right-hand 

rule  dictated  by  the  outward  normal  (integer) . N^ 
must  be  the  smallest. 

* Total  number  of  layers  for  this  element, 
k = Layer  number  (integer) 
t = Layer  thickness 

r = cord  versus  matrix  area  fraction  per  inch 
E^  = Matrix  Young's  Modulus 
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Matrix  Poisson  Ratio 


Ec  - Cord  Young's  Modulus 

V = Cord  Poisson  Ration 
c 

b = Bias  multiplier 
SR  = Rubber  mass  density 
Sc  = Cord  mass  density 

For  identification  purposes  the  letter  E must  precede  the  nodal 
information  and  the  letter  L must  precede  the  layer  flexibility  infor- 
mation. 

5.5  Prescribed  Forces 

At  a given  node  one  may  specify  three  external  force  components. 
The  input  records  read  as 
BEGIN  FORCE  DATA 
k P,  P2  P3 
END  FORCE  DATA 
/ 

/ 

where 

k = Node  number  (integer) 

P^  = Prescribed  force  components  (floating  p.) 

Ip  case  of  pressure  loading  the  normal  (third)  component  is 
assumed  to  follow  the  local  normal  to  the  element  under  consideration. 

5 . 6 Prescribed  Displacements 

At  each  node  the  three  rectilinear  and  two  rotational  displace- 
ment components  may  be  restrained  (i.e.,  equal  to  zero).  These  components 
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are  referred  to  the  base  vectors  of  the  undeformed  reference  surface 
and  labeled  according  to  the  following  schedule,  for  a rotationally  sym- 
metric problem: 

= Component  along  the  parallel 
q2  = Component  along  the  meridian 
q^  = Component  along  the  normal 
q^  = Rotation  component  along  the  parallel 
q^  * Rotation  component  along  the  meridian 

The  input  records  are  exhibited  as  follows: 

BEGIN  DISPLACEMENT  DATA  / 
k M Nx  N2  ...  Njj  / 

END  DISPLACEMENT  DATA  / 

where 

k = Node  number  (integer) 

M = Total  number  of  displacements  restrained  at  the  node 
under  consideration  (integer) 

N^,N2>*..  = The  number  of  displacements  being  restrained  (integer) 
The  contact  problem  is  no  longer  rotationally  symmetric,  and  since 
one  is  dealing  with  a load  of  fixed  direction  the  rectilinear  displace- 
ments are  referred  to  the  base  vectors  of  a fixed  cartesian  frame, 
while  the  rotations  are  decomposed  along  the  base  vectors  of  the  shell 
reference  surface.  Thus, 

q^  = Component  along  the  axis 
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= Cr^ponent  along  the  axis 

q = Component  along  the  x.  axis 

q = Rotation  components  along  the  parallel 
4 

q^  = Rotation  component  along  the  meridian 
The  input  records  are  exhibited  as  follows: 
BEGIN  CARTESIAN  DATA  / 
k M Nl  n2  ...  Njj  / 


k = Node  number  (integer) 

M = Total  number  of  displacements  restrained  at  the  node  under 
consideration  (integer) 

N , N . . . = The  number  of  displacements  being  restrained  (integer) 
1 2 

5.7  Curvefitting  the  Meridian  Reference  Surface 

The  Cartesian  coordinates  of  the  meridian  section  are  defined 
in  this  record.  Thus, 

BEGIN  CURVEFIT  RHOS  / 

X1  x2  ^ 

END  CURVEFIT  RHOS  / 


where 


i = Sequence  number  of  data  points,  (i  = 1,2,...  NPRHO) , (integer). 


x^^2  = Cartesian  coordinates  (floating  p.) 


5.8  Increme 


In  this  record  the  increment  numbers  with  the  corresponding 
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incremental  pressure  are  defined  as  follows: 
BEGIN  INCREMENT  DATA  / 


P1  1 


i P.  / 
1 

™ PIN/ 


where 


END  INCREMENT  DATA  / 


i = Increment  number 


P = Incremental  pressure 


5.9  Contact  Data 

The  candidate  nodes  for  contact  are  arranged  in  a ring-like  fash- 
ion. For  each  ring  one  assigns  the  corresponding  nodes  as  follows: 

BEGIN  CONTACT  DATA  / 

11  1 / 

I Mx  M2  / 

END  CONTACT  DATA  / 


where 


I = Ring  number 

Nj  = Number  of  nodes  in  the  ring 

M , M9>...  Mn  = The  node  numbers  in  the  ring  under  consideration 

1 2 A’  j 


5.10  Print  Options 


The  output  information  is  governed  by  the  following  control  cards: 
BEGIN  PRINT  OPTIONS  / 

ALL  / 


a, 


CONTROL  PARAMETERS  / 

NODAL  DATA  / 

ELEMENT  DATA  / 

CURVEFIT  RHOS  / 

LOCAL  UNIT  VECTORS  / 

NODAL  OUTPUT  TABLE  / 

INCREMENT  DATA  / 

CONTACT  DATA  / 

END  PRINT  OPTIONS  / 

/ 

/ 

END  DATA  INPUT  / 

5.11  Commenting  the  Input  Records 

Following  the  slash  on  the  input  data  cards,  comments  may  be 
inserted.  These  comments  may  be  continued  on  any  number  of  cards, 
having  a slash  for  the  first  character.  Thus, 


BEGIN  NODAL  DATA  / 


1 

0.5 

0.8 

-0.6 

/ 

UPDATED, 

NOV.  11,  1971 

2 

0.6 

1.0 

-0.5 

/ 

BERTRAND 

RUSSELL  DESCRIBED 

/ THE  MATHEMATICIAN  AS  ONE  WHO  NEITHER  KNOWS 
/ WHAT  HE  IS  TALKING  ABOUT  NOR  CARES  WHAT 
/ HE  SAYS  IS  TRUE. 

3 0.8  2.0  0.5  / 


38 


BEGIN  ELEMENT  DATA  / 

1 2 3 10  / THE  ENGINEER  SOMETIMES 

/ PRIDES  HIMSELF  ON  BEING  THE  MAN  WHO  CAN  DO 
/ FOR  A REASONABLE  COST  WHAT  ANOTHER 
/ WOULD  EXPEND  A FORTUNE  ON,  IF  INDEED 
/ HE  COULD  DO  IT  AT  ALL. 

2 4 5 11  / 

• 

• 

5.12  Input  Table  Summary 

The  Input  data  card  set-up  is  summarized  in  this  section 
/ DATA  SET  NUMBER  1 
/ PREPARED  10/1/72 
/ AIRCRAFT  TIRE 
/ 

/ 

BEGIN  DATA  INPUT  / 

BEGIN  CONTROL  PARAMETERS  / 

NODES  N / 

ELEMENTS  E / 

RADIUS  R / 

NRHO  M / 

NPRHO  K / 

LAYERS  L / 

END  CONTROL  PARAMETERS  / 
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. .r.v  a-W.  ^. . . - .......  __ _ ...  . 


/ 


/ 

BEGIN  NODAL  DATA 

i Q <J>  / 

END  NODAL  DATA  / 

/ 

/ 

BEGIN  ELEMENT  DATA  / 

E j Nx  N2  N3  N^  / 

L k t r Er  \ Ec  Vc  b Sc  SR  f / 

END  ELEMENT  DATA  / 

/ 

/ 

BEGIN  FORCE  DATA  / 

k P1  P2  P3  ' 

END  FORCE  DATA  / 

/ 

/ 

BEGIN  DISPLACEMENT  DATA  / 

k M Nx  N2  N3  / 

END  DISPLACEMENT  DATA  / 

/ 

/ 

BEGIN  CARTESIAN  DATA  / 
k M Nx  N2  N3  / 

END  CARTESIAN  DATA  / 
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tr  ■ i'  in  m n 


/ 

/ 

BEGIN  CURVEFIT  RHOS  / 
i xx  x2 

END  CURVEFIT  RHOS  / 

/ 

/ 

BEGIN  CURVEFIT  BETAS  / 

1 xi  x2  ' 

END  CURVEFIT  BETAS  / 

/ 

/ 

BEGIN  INCREMENT  DATA  / 
I P / 

END  INCREMENT  DATA  / 

/ 

/ 

BEGIN  CONTACT  DATA  / 

I N Mx  M2  M3  / 

END  CONTACT  DATA  / 

/ 

/ 

BEGIN  PRINT  OPTIONS  / 
ALL  / 

CONTROL  PARAMETERS  / 
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NODAL  DATA  / 

ELEMENT  DATA  / 

CURVEFIT  RHOS  / 

LOCAL  UNIT  VECTORS  / 

NODAL  OUTPUT  TABLE  / 

ELEMENT  OUTPUT  TABLE  / 

END  PRINT  OPTIONS  / 

/ 

/ 

END  DATA  INPUT 

For  future  extension,  the  code  angle  variation  may  be  described  by 
experimental  data  points.  This  phase  of  the  code  is  not  yet  imple- 
mented, however,  as  the  corresponding  input  data  must  be  present. 
Thus,  augment  the  control  parameter  block  by 
NBETA  2 / 

NPBETA  2 / 
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6.  OUTPUT  DESCRIPTION 


During  the  data  processing  phase,  associated  with  geometrical 
characterization,  the  user  may  exercise  the  print  options  described 
in  Section  5.  The  corresponding  output  information  contains  the  fol- 
lowing records: 

• Input  Data 

The  user’s  input  is  listed.  To  each  card  image  a record 
number  and  a card  sequence  number  are  assigned  for  error 
detection  purposes. 

• Control  Parameters 

0 Curvelinear  Coordinates  of  the  Nodes 

• Element  and  Layer  Data 

• Cartesian  Coordinates  of  the  Data  Points  of  the  Reference 
Meridian 

• Contact  Candidate  Nodes 

• Increment  Data  for  Nonlinear  Inflation 

• Nodal  Output  Table 

• Element  Output  Table 

• Local  Unit  Vectors 

The  ouput  table  of  the  actual  execution  phase  is  not  yet  for- 
malized. Currently,  only  the  generalized  displacements  are  printed 
with  the  element  membrane  stresses  during  the  incremental  inflation 
process  and  tire  rotation.  After  each  step  the  geometry  is  updated 
and  thus  the  Element  Output  Table  and  Local  Unit  Vectors  Table  is 
recalculated.  For  the  contact  problem,  the  nodal  contact  forces  are 
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printed  at  each  itemation  followed  by  the  total  nodal  contact  forces 
and  the  final  geometrical  configuration: 


I 

I ' 


• Intermediate  Nodal  Contact  Forces  carrying  the  title  of 
Solution  Matrix 

• Final  Contact  Forces  carrying  the  title  of  Contact  Forces 

• Final  node  positions  in  a rectangular  cartesian  and  cylin- 
drical coordinate  system 

Thus,  the  code  is  yet  to  be  implemented  by  a complete  output 
module  to  allow  the  analyst  to  select  elements  of  design  interest  for 
stress  and  strain  calculation  purposes. 
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7 . ERROR  EXITS 


Extensive  input  error  checks  are  provided  in  the  data  preparation 
phase.  Each  input  card  is  traced  according  to  its  sequence  number  in 
the  input  deck.  For  cross  reference,  appropriate  record  numbers  (Nodal 
Data,  for  instance  ) are  also  assigned  to  the  input  cards. 

Consider  for  instance  an  erroneous  card  in  the  Element  Data 
(Record  31,  say)  where  identical  node  numbers  are  assigned  to  distinct 

nodes : 

E 1 1 2 2 2 / 

The  error  message  reads: 

RECORD  31)  £ 1 1 2 2 2 / CARD  43 

IN  THE  ABOVE  ELEMENT  CARD,  TWO  OR  MORE  OF  THE 
ELEMENT  NODES  ARE  EQUAL.  ELEMENT  DEFINITION  IGNORED 
Errors  of  this  nature  are  summarized  at  the  end  of  the  data 
processing  phase.  For  the  user's  convenience,  an  input  data  set  is  being 
constructed  which  will  contain  all  possible  logical  errors  with  the  ap- 
propriate  error  diagnostic. 
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n.  Wl-ClaiW., 


_ - - ...  ..  , 


8.  TIMING  AND  STORAGE 


For  large  problems  it  is  important  to  estimate  the  needed 
execution  time  for  both  central  processor  and  peripherial  operations. 

At  this  time  not  enough  data  is  available  to  either  construct  appro- 
priate formulas  in  terms  of  major  computational  parameters  or  graphs 
based  on  direct  experiemtation.  On  the  CDC  6600  machine  at  the  WPAFB 
under  the  RUN  system  the  central  processor  time  may  be  estimated 
according  to  the  formula 

CP  (minutes)  =10  3 * (E  * L * S) 

where 

E = Number  of  elements 

L = Number  of  layers 

S = Number  of  steps  for  incremental  inflation  . 

For  large  problems  the  peripherial  time  is  roughly  that  re- 
quired by  the  central  processor. 

The  program  does  not  yet  contain  output  information  for  min 
imum  execution  field  length  requirements  during  the  loading  phase.  If 
the  declared  field  length  is  too  small  for  execution,  an  allocation 
error  message  will  appear  at  the  corresponding  phase  of  the  code.  Be- 
cause of  the  size  of  the  fourth  overlay,  substantial  storage  is  re- 
quired even  for  small  (50  elements)  problems,  such  as  120,000  central 
memory  in  octals.  The  largest  test  case  (200  elements)  required 
135,000  central  memory  in  octals. 
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9 . SAMPLE  INPUT 


In  this  section  a pathological  example  is  considered  to 
demonstrate  the  structure  of  the  input  data  which  covers  all  phases 
of  the  computer  code. 

The  problem  under  consideration  is  the  inflation,  rotation, 
and  contact  of  a strip  along  the  meridian  of  a toroidal  shell  shown 
by  Figure  1.  It  is  assumed  that  the  strip  is  of  uniform  thickness 
consisting  of  two  layers.  The  cord  angle  varies  along  the  meridian 
according  to  the  classical  lift  equation  [1].  The  input  set  listed 
below  is  annotated  for  clarity  in  presentation. 
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NONLINlAR  STRlp 


/ 

/ 

/ 

/ 

/ 

/ 

crtL.CK.Oo  I / TuKNS  ON  PRINTS  cP"  RcollNco  NmMcC 


/ 

/ 

BEGIN  DATA  INPUT  / 

BEGIN  CONTROL  PARAMETERS 
ELEMENTS  12  / 

NODES  11  / 

RAD  I Ub  9.15  / 

NR HO  10  / 

NPKHO  27  / 


NPETA  2 / 

NPBETA  2 / 

LAYERS  2 / 

GREEN  0.95 
SPEED  100. 
NONLINEAR  / 
I NCREMENT  b 


NOT  ACTIVE 
NOT  ACTIVE 

/ 

/ 

5 / 


/ NO  PARTICULAR  ORDER  IS  ASSIGNED 


BUT  MUST  BE  PRESENT 
BUT  MUST  BE  PRESENT 


CONTACT  -0.05  vouu.  / WHEEL  LOAD  IS  NOT  ACTIVE 

NORING  3 / 

NO DM AX  3 / 

end  control  parameters  / 

BEGIN  NODAL  DATA  / 

1 o.  0.  / GIVEN  IN  RADIANS 

2 0.^2  0.  / 

3 0.01  0.056  / 


A 0.  0.112  / 

5 0.02  0.112  / 

6 0.01  0.16b  / 

7 0.  0.224  / 

6 0.02  0 , 22  A / 

9 0 . 0 1 o,28  / 

10  o.  0.336  / 

11  0.02  0. 3 J6  / 

END  NODAL  DATA  / 

/ 

/ 

BEGIN  ELEMENT  DATA  / 


E 1 1 2 3 2 / 


L 

1 

0 .0645  0.435 

450* 

0*49 

150c00* 

0*  7 

i 

• u • Qo  0 1 0# juol 

o 

• 

/ 

L 

2 

0 • 06 A5  0.435 

45  0* 

0*  A9 

1 5ooOO • 

0*7 

“ 

1 • u • \j  0 0 1 u*  uJul 

0. 

/ 

E 

2 2 5 3 

2 / 

L 

1 

0*0643  0.433 

450  • 

0*Av 

15ocoG* 

0*  7 

i 

• 0«0001  U«UUol 

0 

• 

/ 

L 

2 

0*0645  0.455 

450  • 

0*^Y 

1 5 b 0 0 U . 

0*7 

- 

1 • 0 • 0 00 1 0 • U U 0 1 

0. 

/ 

E 

3 1 3 A 2 

/ 

L 

1 

0 • 0645  0.455 

450* 

0*49 

156000* 

0*7 

1 

• o.oooi  o.ouui 

0 

• 

/ 

L 

2 

0*0645  0.455 

450  • 

0*49 

1 56ou0 • 

0*7 

“ 

1*  uQu  1 Ui  c/OO  1 

0* 

/ 

E 

4 3 5 4 

2 / 

L 

i 

0.00*43  0.405 

*4  0 0 • 

0*4  3 

c.7 

1 

• OwU  1 U*uv/Ul 

0 

• 

/ 

L 

2 

0 • 0 o 45  0 . 4o 5 

430  • 

0*^5 

1 5buOO • 

0.7 

- 

1#  u • 00 0 1 u • o U j 1 

0. 

/ 

E 

5 4 5 6 

2 / 

L 1 0*0643  0.43b  4bO»  0*44 
L 2 0*0643  0.435  430*  0*44 
E 6 b 8 6 2 / 

L 1 0 • 0 64b  0 • 4 3 3 4 30  • 0*49 
L 2 0 • 0 o4  3 0.43  b 430.  0«44 
E 7 4 6 7 2 / 

L 1 0*  064b  0 . 4 3 b 4bO.  0»4^ 
L 2 0 • 064b  0.4 Jb  430*  0 • 4 4 
E 8 6 8 7 2 / 

L 1 0.064b  0 • 4 j 3 430*  0*44 
L 2 0.0  64  b 0 . 4 b b 4b0.  0.44 
£9  7 8 9 2 / 

L 1 0 » 0 1>4  3 U.4^3  4 3 0 • 0*4/ 

L 4 0.0o43  0.433  430.  0.4? 

E 10  8 11  9 2 / 

L 1 0 • 0 0 4 3 0.433  430*  0.4* 
L 2 0*0643  0.433  430*  0.49 
E 11  7 9 10  2 / 

L 1 0*0643  0-433  450*  0*49 
L 2 0*0643  0.433  430*  0*49 
E 12  9 11  lo  2 / 

L 1 0*0643  0.43b  4bO*  0*44 

L 2 0*0343  0.433  430*  0*44 

END  ELEMENT  DaTm  / 

/ 

/ 


BEGIN 

FORCE 

DATA 

/ 

1 

0. 

0. 

1. 

/ 

NOT 

ACTIVE  BUT  MUST  BE  PRESENT 

2 

u. 

0. 

1. 

/ 

SEE 

PRESF  SUBROUTINE  FOR  CONCENTRATED  LOADS 

3 

0. 

0. 

1. 

/ 

4 

0. 

0. 

1 

X • 

/ 

5 

0. 

0. 

1. 

/ 

6 

0. 

0. 

1 . 

/ 

7 

0. 

0. 

1 . 

/ 

8 

0. 

0. 

1. 

/ 

9 

0. 

0. 

1 . 

/ 

10 

3 , 

0 

. 1 . 

/ 

11 

u . 

0 

. 1 . 

/ 

END  FOKCE  DATA  / 


/ 

/ 


f 

BEG  I N 1 

D13PLAC 

EMcNT  JAT A / 

1 

5 

1 2 3 

4 3 

/ 

FIXED 

2 

3 

1 2 3 

4 3 

/ 

FIXED 

3 

2 

1 

3 

/ ZERO  PARALLEL 

DISPLACEMENT 

AND 

MERIDIAN  ROTATION 

4 

2 

1 

b 

/ 

5 

2 

1 

3 

/ 

6 

2 

1 

5 

/ 

7 

2 

1 

5 

/ 

8 

2 

1 

3 

/ 

9 

2 

1 

5 

/ 

10 

b 

1 2 3 

4 3 

/ 

11 

5 

1 2 3 

4 5 

/ 

END  DIoPLACEMtN I DATA  / 
/ 

/ 


1 3 oOb U « 0*7 

1.  U.0001  O.Ooul  ( 

0. 

/ 

1 56  J JO  » 0*7 

“1.  0.0001  J • j 0 0 1 

0. 

/ 

1 36 J^o » j»7 

1*  o • u J U 1 0 . U u u 1 i 

o. 

/ 

1 36 uu u » U»7 

— I.  o • oOOl  u • o u U 1 

0. 

/ 

1 5oouO  . J * 7 

1.  u.uoOl  u.uoul 

0 • 

/ 

1 36UUU  * 0*7 

-1.  O.UOOl  O.jOOI 

o. 

/ 

1 boUuO  » 0*7 

1.  0.0001  0.0001 

o. 

/ 

156 jjO  * U»  7 

“1.  o.uOOl  O.ejOI 

o. 

/ 

lsoU-ju.  \J  . / 

1.  U.OuOl  u.u^ul 

0. 

/ 

i3u juo.  0»  7 

~i.  u.oOOl  O.uOOl 

0. 

/ 

lbooou.  0*7 

1.  u.OoOl  O.JoUl 

0. 

/ 

136UU0.  0*7 

-1.  o.oool  0.0001 

0. 

/ 

1 56Cl>0  • 0.7 

1.  0. OoOl -0.0001 

o. 

/ 

136 OUU*  0*7 

— 1*  u • u 00 1 j.ubOl 

o. 

/ 

1 3b ju u • 0*7 

1 . u. OuOl  O.uuul 

0. 

/ 

1 b 6 ou  0 » U»7 

— I*  ^ . 0 0 u 1 u « o J 0 1 

o. 

/ 
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BEGIN  CARTESIAN  DATA  / 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


5 

3 

2 

2 

2 

v 

2 

2 

2 

5 

3 


1 

1 

2 

2 

2 

2 

2 

2 

2 

1 

1 


2 

2 

5 

3 

5 

3 

5 

3 

5 

2 

2 


END  CAR  T E 4 I AN  DATA  / 
/ 

/ 

bcGIN  CuRVEFlT  RH03 


1 

4.69  0.  / 

2 

4.68  0.37 

/ 

3 

4.66  0.74 

/ 

4 

4.63  1.1  / 

3 

4. do  1.47 

/ 

6 

4.51  1.94 

/ 

7 

4.4  2.19  / 

8 

4.25  A. 32 

/ 

9 

4.03  2.84 

/ 

10 

3.63  3.13 

/ 

11 

3.57  3.39 

/ 

12 

3.29  3.62 

/ 

13 

2.97  3.81 

/ 

14 

2.63  3. 9 '3 

/ 

13 

2.29  4.03 

/ 

16 

1.89  4.04 

/ 

17 

1.52  3.97 

/ 

18 

1.16  3.8o 

/ 

19 

0.82  3.73 

/ 

20 

0.5  3.55 

/ 

21 

v . lo  3 .33 

/ 

2 2 

o.  3.24  / 

2 3 

-0.25  3.1 

/ 

24 

-0.425  3.C 

) / 

25 

-0.6  2.9 

/ 

26 

-0.75  2.6 

/ 

27 

-0.95  2.7 

/ 

ND 

CURVFIT  RHQj 

/ 

dEGIN  increment 

DATA 

1 

2.  / 

2 

2.  / 

3 

2.  / 

4 

10.  / 

5 

20.  / 

3 

3 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

3 

3 


45/  FIXED 
43/  FIXED 
ZERO  Y-DISPLACEMENT/MERIDIAN  ROTATION 


4 

4 


FICTICIOUS  POINTS  TO  AVOID  OSCILLATION 


END  INCREMENT  i) ATA  / 


BEGIN  CURVFIT  BtTAS  / 
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1 1.  2.  / 

2 1.  2.  / 

LND  CUkVUT  ocTmS 


NOT  ACTIVE,  BUT  MUST  BE  PRESENT 
NOT  ACTIVE,  BUT  MUST  BE  PRESENT 


7 


BEGIN  CON  F ACT  D a T A / 
111/  THE 

2 3 2 3 A / 

3 2 i>  6 / 

END  CONTACT  DAIa  / 

/ 

/ 

bEGIN  PRINT  OPTIONS  / 
CONTROL  PARAMETERS  / 
NODAL  DATA  / 
element  DATA  / 

CORVFIT  RrlOS  / 

CURVFIT  oE  t as  / 

-GCAL  UNIT  VECTukS  / 
NODAL  OUTPUT  TAuLE  / 
INCREMENT  DATA  / 

ELEMENT  OUTPUT  TABLE  / 
END  PRINT  OPTIONS  / 

/ 

/ 

END  DA  f A INPUT  / 


FIRST  CONTACT  NODE 


MUST  BE  LABELLED 


1 
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10.  PROGRAM  LISTING 


The  Load  Map  presented  below  was  generated  on  a CDC-6600 
machine  using  the  RUN  compiler  at  the  Wright-Patterson  Air  Force 

I 

Base  installation.  The  corresponding  program  listing  is  avail- 

I 

able  upon  request  for  qualified  applicants. 

I I 

i 


bo 

bo 


fa* 

w 

fa* 


fa* 

fa* 


I 

o I 
o h- 
o <J 
o Z 
O UJ 

O _J 
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z 
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o o 
O O 
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^ O 


ro 
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in 
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£ o 
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UJ 

o z 
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CL 
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a z o' 
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CD 
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►•aw 

Z > -H 
O h (M 
O • t-1 


> UJ 
< Q 
-J  O 
O'  £ 
UJ 

> a 

o < 


i 

in 

M/IN 

03  UJ  -rH 

n-  or  -» 

O Q rl 
W O o 

H < O 


N.4-N*JHh.Sin 

<-vD*Hroin«n*-ioj 

oHWwwojroin 

inininminininin 

t\JMCMCM(MC\J(M(VJ 

a o o o o o o o 


■Jinr((\Jr(NOOrO^KNMMJ3r(rOtHON 
J JMNJinO  JiflOd  NN-ONWrONN  f~*  h- 

lAiD^NKooJiniMnsriHNroj  o in  k 

ininininino>ijD(j3«nkOvOvON-N«ooooo^ 
{M(M(MC\J(NJ(MN(M(M(MnjCM(M<\irororOrOIOIO 
o o o ooooooooooaooao  o o 


i 


t£  I 
UJ  I 
a c 


N I 

* 


c,  a 
-i  c 
a uj 
< t*  c* 

X C.  M 

u.  I h- 


haihoauz  h 
aDQ-orwow 
awrouozaw 

h-KS-h-h-k-h-  Qi 

<<<J<<<1<Q ' 

rrir rr r w 


in 


r i/io  »-  jfflh 
3 Q O < 3 D < 

1/1  Of  < r £ (/)  U WUI  ru.rnu.BiuU4W 

aaOoooojhNOHO(flujhuj»-Q  o 

aZUJlJWUJJlLl0l/)Q.ZO>>D>U)DH 


o 

O' 

O UJ 


>-  o 

r < o o 

_ UJ  -I  h-  -I  < O' 

hlLhaaKfflklM 


az  lii  UJUJUJ  JU.W  WQ.  ZQ>>3>IUQH 

uj  h > >>><^r  won  jwoooo^ifl 


a 

< 

a: 


o 

o 


54 





UK  * - ' " •'•w-  ■ 


N. 

J- 

ro 

OJ 

0 

lO 

OJ 

OJ 

o- 

ro 

OJ 

vO 

OJ  LO 

0 

OJ 

OJ 

«H 

«H 

N- 

0 

ro 

0 

ro 

J-  -T 

0 

J' 

0 

ro 

J- 

LO 

05 

O 

a? 

05 

•H  «H 

lD 

^4 

o- 

0 

O 

H 

«H 

LO 

OJ 

OJ 

OJ 

OJ 

O O 

OJ 

O 

OJ 

ro 

ro 

ro 

O 

OJ 

0 

0 

0 

O 

O O 

0 

O 

0 

0 

0 

0 

O 

O 

>• 

to  «t  o 
z jhq; 
►-  or  a ljj 

H liJK  3 

z > =>  o 

►H  O O ^ 


X 
lli  a 
ar  »- 

h*  <1 

►-  z 


i 

to 

to 


_ JLflKLf'NOOiDLf»iDf\jrO(\JO*HHlf\HO 

UJCJLOO^rOLr'N.rjr'JJ-<'r>-r\JLr'»H-J’LrirjPJO 


w _ . . . ..  < J-  ITi  r\J  (M  <x> 

C£-J’O*-40J(\|CJ0Jf0L0Lf\LCLDN*r,»OOOWCJ.J' 

Dt-»LnLnLnLALiMALiM^Lr'LA^Ln^^^^^^^ 

□ 0(\jt\Jt\JfMN(\JNN(\Jf\jry(V  (\J  OJ  (\J  CJ  OJ  OJ  OJ 

<0000000000  000  0000000 


<0 

ou 


LO 

lD 

OJ 


LO 

Jf 

S J3  S r( 

ro  0 

0 ro 

O 

■J  n n 

N. 

IO 

CJ  J-  ^5  O 

LO  OJ 

OJ  OJ 

ro 

O-  OJ  *4 

LO 

05 

LO  LO  N-  O 

0 

»*♦ 

ro 

OJ  ro  OJ 

05 

05 

kC  iD  iD  N 

N»  N» 

fv  N. 

N. 

K S N 

OJ 

OJ 

OJ  OJ  OJ  OJ 

OJ  OJ 

OJ  OJ 

OJ 

OJ  OJ  OJ 

O 

O 

0000 

0 0 

O O 

O 

000 

1 

V 

a 

r- 

a d 

>- 

0 0-0 

z 

to 

0 ►“ 

_J 

CD 

UJ  UJ 

►-  O'  to  to 

0 

<r 

to 

X 

V 

Z 

0 a 
z z 

4- 

O' 

CL 

0 Z UJ 

< 

LlI  Z> 

0 

0 

< 

Z> 

3 

000 

O' 

OJ 

ro 

•j  tr 

z 

0 

U_  Ckr 

UJ 

UJ 

UJ 

Z 

UJ 

CL 

to 

Z 

< to  z 

a 

to  co 

O' 

<r 

z 

z 

to 

<_> 

O CL  CL 

0 

UJ 

0 

0 

0 cr* 

z 

Z L- 

O ►- 

H* 

u- 

K5* 

tu 

1 

cxl 

U* 

H- 

►- 

H- 

ct  ar  a 

0 

CJ 

0 

CJ 

0 

z u.  t* 

-J 

►“ 

N» 

0 

0 

000 

UJ 

»-t  O 

0 z 

CO 

to 

to 

0 

M 

<r 

<t 

<r  < <3 

< 

CkT  CkT 

z 

UJ 

UJ 

UJ 

UJ 

_J 

CL  UJ  -J 

u. 

to 

to 

CL 

CL 

CL 

a cl 

a. 

r-  O 

-J  <50 

>• 

>• 

>-  £ 

1 

U* 

Z 

Z 

z 

z z z z 

LU  UJ 

►H 

> 

> 

> 

> 

OOO 

V 

z 

to 

O 

O 

O 

000 

to  -J 

X Cf 

to 

to 

to 

LU 

55 


J 


MATSU8  025206 
MATRAN  025344 
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